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ABSTRACT 



Least squares methods are used with NIMBUS IV SIRS-B 
clear-column radiances as independent variables and two 
identical-period sources of geopotential heights (at the 
500 mb and 300 mb pressure levels) as dependent variables. 
Regression equations are developed for each of the sets of 
geopotential heights for three latitude bands in the 
Northern Hemisphere. Three-day and four -day data bases were 
used on each regression specification. Each set of equa- 
tions was tested on independent data at a subsequent 
composite period of 12-hours and 24-hours following the 
dependent pooled samples. The regression results of both 
the dependent and independent tests were examined for 
possible operational usefulness. 
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I. INTRODUCTION 



The Satellite Infrared Spectrometer (SIRS-B) carried on 
the NIMBUS IV spacecraft system relays to National 
Aeronautics and Space Administration (NASA) , at Goddard 
Space Flight Center, radiance observations taken in seven 
distinct channels of the 15 ym band of carbon dioxide, and 
one at 899.0 cra"^ in the atmospheric window. In addition, 
there are six channels in the water-vapor rotational band 
also being relayed by NIMBUS IV, but these latter six 
channels will not bear on the subject of this dissertation. 

Each of the eight primary channels (699.3 cm - ^, 

677.8 cm 692.3 cm \ 699.3 cm"'* - , 706.3 cm \ 714.3 cm -- * - , 
750.0 cm \ and 899.3 cm '*’) is sensitive to temperature 
variations in the tropospheric-stratospheric range of heights. 
In addition, each channel has a maximum transmissivity in a 
discrete height range above the surface. Smith, Woolf, and 
Jacob (1970) presented a graph depicting dT /dlnp versus 
pressure where t is the transmissivity and the derivative, 
dT /dlnp, becomes a weight factor of the relative contribu- 
tions of the various elevations to the total emitted 
radiance in the channel. It can be seen in Figure 1 and 
Equation (1) , after Smith et al (1970) , that the different 
channel radiances are most responsive to changes in tempera- 
ture which occur at the level of their peaks in dT /din p. 
Each channel tends to signal the occurrence of temperature 
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changes at or near its own pressure-level peak, as 
discernible from Figure 1. 




Fig. 1 Derivative of transmittance with respect 
to the logarithm of pressure for eight 
SIRS-A channels [after Smith et al 
( 1970 )] . 



The radiative transfer integral applied to a clear 
column of air for which the sounding T(p) is known gives 
the radiance as 



N 



~ ( Ps> 

(Vi) = B^Vi»T(p s )J T( Vi,p s ) - j B [Vi-T(p)J dTi 



( 1 ) 



where B 



[Vi.T(p)J 



is the Planck radiance at wave number 



Vi" and T(p), the temperature at pressure level p. In 
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indicates a surface value and 



Equation (1) the subscript "s" 
t(Vl» P) is the fractional transmittance of the atmosphere 
in the spectral interval centered at from pressure 

level p up to the satellite. If the sounding is known, 
the radiance in the channel centered at wave number 
can be computed since t( V^, p) is known as a function of 
pressure p Smith (1970) . The inverse problem has become 

more relevant since the advent of the SIRS instrument, 
namely in the use of observations in the eight channels to 
deduce the clear atmosphere structure T = T(p) . The usual 
method of solving an integral equation for the function 
B( Vi# T) was found to be unstable, but it was found £ Smith 
et al (1970)J that regression methods relating temperatures 
and geopotential heights at standard pressure levels were 
quite well described by the results of SIRS scan-soundings 

r 

Smith, Woolf, and Jacob (1970) found it convenient to 
convert the clear column radiances into equivalent 
"brightness temperatures," Tg ( Vi). The definition of 
Tb( ) is based upon the particular Planck function for 
temperature t b( Vi) which equals the clear-column channel 
radiance N( Vi) , that is 

Epi- t E (Vi’] - N(Vj.> ( 2 ) 

so that by solving for T B we have 

T B (Vi) = CsV^i | In £ C 1 V i 3 /N( V ± )J + 1 | 1 (3) 
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Here the constants and C 2 have the values 

= 1.19061 x 10 erg-cm -sec A -steradian 
C 2 = 1.43868 cm-°K 

Smith et al (1970) derived regression equations for contour 
heights at standard pressure levels (see Table 1) utilizing 
the eight simultaneous brightness temperatures Tg(VjJ as 
independent variables. As dependent variables they employed 
a large sample of radiosonde-reported geopotential-heights 
at the mandatory levels. 'Their results for the month of 
September, 1969, using a two-week sample, was stated in the 

nonlinear regression form 

8 

Z( Pj ) = Z(Pj) +' X/ b ( V'i'Pj).[TB(-^i) ~ Tgtv'i)] + 

1=1 (4 

8 

X>’(wPj) [Wi) - Wi )] 2 

i = 1 L 

with an analogous regression equation for T(pj) . Their 
regression coefficients for the September, 1969 SIRS-A 
example are shown in Table 1, where the coefficients for 
T (pj ) are also listed. For the data period, 24-29 December 
1970, employed in this study, an analogous set of regression 
equations based upon updated SIRS -climatology would be used. 

It is clear that the form in Equation (4) , as employed 
by Smith et al (1970) , essentially builds from the tempera- 
ture structure of the remotely -sensed atmospheric columns 
v^ith p s = 1000 mb and reconstructs a height profile by 
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ecsxe 



statistically relating Z(pj) through the hydrostatic 
equation to Z(pj.) at any other level desired. 

The objective of this study was to obtain statistical 
relationships for geopotential heights of the 500 mb and 
300 mb pressure levels by stepwise regression analysis 
between these variables and the eight clear-column satellite 
radiances. Further, it was desired to do this without prior 
computation of the brightness temperatures. The latter 
procedure is an unnecessary step if corrected radiances are 
known, as was the case in this study. If a sixteen- 
predictor equation were developed. Equation (4) would then 
be modified to include the N^ and N^ terms (i = 1,...8) 
with appropriate coefficients. 

With the use of NIMBUS IV SIRS-B corrected radiances, 

Nj[, and National Environmental Satellite Center (NESC) 
regression-computed contour heights at the same geographic 
locations, diagnostic regression equations for the 500 mb 
and 300 mb geopotential heights were developed. Separate 
regression equations were developed using both three-day 
and four-day pooled data bases. The regression equations 
developed were tested for validity on an independent sample 
in each case. 

The same procedure was applied to Fleet Numerical 
Weather Central (FNWC) analyzed geopotential heights at the 
500 mb and 300 mb levels for the same positions and times. 

In the latter case, the contour heights had in no way been 

related to SIRS-B scan-soundings. Thus, any diagnostic and 
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prognostic relationships discovered here would have to be 
based upon the statistical correlation coefficients 
R [Z(5, 3) /NjJ , i = 1,...,8, as previously suggested by 
Smith et al. 



18 



II. DATA PROCESSING 



A. THE ORIGINAL NESC DATA 

The original data obtained from the National Environmental 
Satellite Center (NESC) consisted of a computer listing, by 
days, of selected SIRS-B scan spots during the period 
Dec. 1970 - Jan. 1971. Included in this listing were the 
geographical sites and times of each observation, amounts 
and heights of effective cloud cover, temperatures and 
heights of the standard levels, and a double set of 
radiances. The latter were the original set of uncorrected 
radiances and the cloud-corrected values for each channel.- 
The uncorrected radiances in the window channel were missing 
since that channel had become inoperative by December, 1970. 
The sample of scan-spots presented by NESC was therefore 
restricted to sites having Known surface temperatures so 
that a realistic effective window channel radiance, B(v",T s ) 
for a surface pressure level p s = 1000 mb could be assigned 
for each location. 

It should be pointed out at this stage that the SIRS-B 
radiometer carried aboard the NIMBUS IV differed from the 
SIRS-A radiometer aboard the NIMBUS III in that the SIRS-B 
radiometer had a "side -looking" capability as well as the 
vertical mode. This permitted a wider range of so landings 
in SIRS-B than was possible with SIRS-A. These side-view 
radiances v/ere normalized to a single atmospheric 
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emitter-depth located over the scan spot so that all the 
data was compatible, apart from variations in cloud cover 
and surface-layer lapse variations. 

Because of the scan-spot area of (225 km)^, most spots 
viewed were partially cloud-filled with some effective 
fraction, A. This fact complicated the determination of the 



cloud -corrected 


T (p) 


profiles since the radiances were 


reduced to 






N(Vi) = A 


| B [ 


Vi' T i<Pc>] T<V i .P c ) 






) 


— 


/ B 
1 


[ VVT(p) J d Ti > 


-r 


(1 - A) 


| B [Vi,T(p s )] j (Vi<P s ) 



- f B [ V'i'Tte)] d T i| 

assuming a one-level cloud model with the top at p c . Here 
the part of Equation (5) within the second brace is that 
part of N(v'i) which represents the desired uncontaminated 
radiance in channel i which escapes through the cloud-free 
fraction, 1 - A. 

Smith et al (1970) have outlined a method for determining 
a cloud-contamination correction which yields the estimate 
for corrected radiance N c (Vi) in each channel. Since the 
work in this study was dependent upon the accuracy of the 

corrected radiances, a brief explanation of this correction 
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procedure is given, based on the one-level cloud model with 
top at level p c . Equation (5) can be broken into two parts, 

N c (Vi) = N ( V ± ) + C( %/.) (6) 

where the correction C (V is designed to eliminate the 
effects of cloud-reduction of the radiance, and may be 
obtained by dividing Equation (5) by the factor 1-A. How- 
ever, this will not be a permissible procedure since the 
cloud cover. A, may be very close to unity. Instead, Smith 
et al (1970) suggest an iterative procedure whereby C ( V 
of Equation (6) may be obtained in the form 

C( V ± ) = A | N c (Vj.) " X |y r P c ,T(p)] | (7) 

The function X ^ v /-l'P c ' T (p) J represents the contents of the 
second brace of Equation (5) . Successive corrections 
C(v^) are then determined for each channel 4,...., 8, by 
making the assumption that the radiances observed in 
channels 2 and 3 (7 50.0 cm"'*' and 714.0 cm"* - ) are correct 
during iterations while correcting the remaining set of 
radiances N^,...Ng. 

In determining a first guess to a clear column atmosphere 
T(p) the five most opaque channels in Figure 1 are assumed 
correct at the first iteration, as well as the pre- 
determined window channel. The NESC correction scheme then 

makes use of a six-predictor regression equation for T(p) 

♦ 

(analogous to Equation ( 4 ) but lacking channel 2 and 3 terms) 
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which then specifies the required first guess T(pj) at all 
standard levels. With this profile now specified the 
best-fit consistent choices of A, 1-A, and p c which comes 
closest to providing the best-fit solutions N( v/^) to 
Equation (5) for channels 2 and 3 are obtained. In this 
test p c is subject to variation over the standard levels 
p = 200, 250 ,....850 mb and the best-fit A, 1-A, are 
thereby chosen. At each iteration, there is an improved set 
of six corrected radiances, which leads to an improved 
specification T(pj) . This then leads to improved newly 
corrected radiances N c ( \A) (i = 4,....,8) in accord with 
Equations (6) and (7) . Finally, in a limited number of 
iterations C(v'j^) approaches a convergent value. At this 
final stage, the six-predictor regression equation for 
T(pj) gives estimates of corrected radiances in channels 2 
and 3. The listed statistics relating p c and A are of 
no particular importance since these quantities relate only 
to an effective cloud cover over an area (225 km)^, and the 
cloud-depth is considered irrelevant in obtaining N c Vi>- 
The important factor to be tested in this paper is the 
concept that geopotential heights as taken from a large 
radiosonde sample are highly correlated with the set of 
clear-column radiances, i.e., to the estimated T(p) 
structure of the atmosphere. Thus the NESC SIRS-B data 
obtained by the December 1970 version of the NESC regression 
equation system provided also the geopotential heights and 



22 



temperatures at all standard levels in the atmosphere from 
p = 1000 mb to 10 mb. 

In developing the NESC regression equations. Smith et al 
(1970) utilized a data sample of approximately two weeks, 
continually updating the sample period closer to the time 
when the independent calculations were made. The heights 
used in their regression equations were interpolated in 
space and time from National Meteorological Center analyzed 
fields at standard pressure levels to the geographical loca- 
tions and observed times of the SIRS measurements. Only 
portions of analyses and NIMBUS III data close to radiosonde 
observations were used for developing their dependent data 
regression coefficients giving heights from corrected 
radiances. These NESC dependent samples were divided into 
latitude bands (18-40°N, 35-55°N, and 50-80°N) . A similar 
technique was used by NESC in making available the data from 
December 1970, a portion of which was used in this thesis. 
This data consisted of corrected radiances at scan-spots 
which were used in conjunction with an eight-predictor 
regression equation for generating temperature and height 
values for the standard levels from the SIRS observations. 
These computed heights and temperatures appeared on the NESC 
data listings for 24-28 December 1970, which were used in 
this study. 

The NESC regression-derived heights of the 500 mb and 
300 mb levels, position, time, and clear-column radiances 
for each sounding from 24-27 December, 1970, were used for 

the regression analysis described in Chapter III. 
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B . FNWC DATA 



Height fields for 0000 GMT and 1200 GMT covering the 
period 24-29 December 1970, were obtained from the data- 
tapes of Fleet Numerical Weather Central (FNWC) at Monterey 
for both the 500 mb and 300 mb levels. The contour heights 
were interpolated in space and time to the same positions 
and times of the satellite data obtained from NESC. The 
spatial values of geopotential height were computed by use 
of a double-quadratic Bessel interpolation scheme applied to 
both the synoptic times preceding and following each satel- 
lite observation. Then the heights were interpolated 
linearly to the time of the satellite observation. The 
computer program used to accomplish the interpolation scheme 
is listed in the Appendix. 

The FNWC data provided an independent set of geopotential 
heights for the 500 mb and 300 mb levels from which a second 
set of regression relations could be derived for comparison 
with the regression analyses performed on the NESC data. 
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III. DEPENDENT DATA REGRESSIONS 



A. DEVELOPMENT OF REGRESSION EQUATIONS USING NESC GENERATED 
VALUES OF 500 MB AND 300 MB GEOPOTENTIAL HEIGHTS. 

The NESC-generated 300 mb and 500 mb pressure levels 
(Z n ( 5) and Z N (3) respectively) were selected for investi- 
gation with the corresponding 500 mb and 300 mb FNWC 
analyses, which presently have no SIRS input. A multiple 
linear regression analysis technique was used in this study 
corresponding to the 8-predictor scheme of Smith, Woolf, and 
Jacob (1970) , but containing no nonlinear terms. The speci- 
fic program used was BIMED 02R, a version of the stepwise 
linear regression technique currently on file at the 
Statistical Library of the W. R. Church Computer Facility 
at the Naval Postgraduate School, after Dixon (1966) . The 
result of this technique when applied to simultaneous data 
samples of the form [Z^/N^,^, . . ,Ng] is to produce 
regression equations of the form 




Here the subscripts 1,..,8 indicate the channel numbering 
system and are also used for variable numbering. A series 
of dependent test regressions we re performed on both Z^(5) 
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and Z n ( 3) using a pooled three -day data base and a pooled 
four-day data base of simultaneous values of SIRS-B correc- 
ted radiances. Due to the nature of the data available, 
each data-day presented a composite of scan-soundings 
centered about 1200 GMT, varying up to six hours on either 
side of this central time. The data was further divided 
into three overlapping latitude bands. The limits of each 
band are listed below: 



Band 


1 . 


20°N 


- 40°N 


Band 


2 


35°N 


- 55°N 


Band 


3 


50°N 


- 70°N 



The resulting coefficients (C 0 , C^,....,Cg) for the three- 
day and four -day regression equations are listed in Tables 2 
and 3 for the 500 mb and 300 mb dependent tests, 
respectively. 



26 



TABLE 2. 500 mb NESC Dependent Regression 

Coefficients for Three- and Four-Day 
Data Base, Giving Best-fit for Eq. (8) 



3-Day Sample Base 





Band 1 


Band 2 


Band 3 


Sample Size 


171 


190 


177 


c o 


5473.77 


5689.80 


4620.48 


C 1 


6 , 66 


-2.80 


-2.90 


c 2 


-17.56 


-3.32 


-11.68 


c 3 


23.85 


41.83 


62.59 


c 4 


28.06 


-10.83 


-17 . 87 


c 5 


-62.22 


- 66.78 


-36.90 


C 6 


-6.80 


4.02 


-4.14 


C 7 


21.20 


19.62 


11.74 


C 8 


-0.33 


-1.40 


* 




4 -Day Sample 


Base 






Band 1 


Band 2 


Band 3 


Sample Size 


171 


227 


216 


c o 


5497.55 


5756.70 


4555.32 


C 1 


6.53 


-1. 18 


-4.45 


c 2 


-15.98 


-7.7 5 


-6.38 


c 3 


20„ 13 


45.24 


54.55 


C 4 


34. 33 


-12.75 


-10.93 


C 5 


- 67.57 


- 64.77 


-35.75 


C 6 


-5.88 


2.64 


-7.94 


C 7 


24. 18 


19.54 


1.00 


c 8 


-3.26 


-1.13 


1.45 



Indicates an insignificant contribution to the analysis 

of variance by this predictor. 
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TABLE 3 



. 300 mb NESC Dependent Regression Coefficients 

for Three* and Four -Day Data Base, Giving 
Best Fit for Eq. (8) 



3-Day Sample Base 





Band 1 


Band 2 


Band 3 


Sample Size 


171 


190 


177 


c o 


8976.59 


8490.63 


7166.52 


C 1 


12.79 


4.00 


1.57 


c 2 


-22.36 


-6.08 


-15.54 


C 3 


4.04 


21.62 


54.11 


C 4 


87.54 


47.07 


23.71 


C 5 


-105.20 


-96. 10 


-48.47 


C 6 


-15.77 


4.64 


-4.34 


C 7 


23.90 


15.88 


6.14 


c 8 


5. 67 


1.62 


1.29 




4 -Day Sample 


Base 




. 


Band 1 


Band 2 


Band 3 


Sample Size 


206 


227 


216 


Co 


8986.32 


8613.84 


7102.40 


C 1 


13. 26 


6.65 


* 


c 2 


-21.76 


-13.75 


-10. 18 


C 3 


* 


29.09 


45.98 


c 4 


95.49 


41.53 


30.92 


c 5 


-110.81 


-93.15 


-47.49 


c 6 


-14.34 


3.13 


-8.48 


C 7 


28.03 


16.26 


5.37 


c 8 


1.25 


1.65 


2.94 



★ 



Indicates an insignificant contribution to the analysis 
of variance by this predictor. 
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The coefficients in Tables 2 and 3 are considerably 
different in each latitude band which tends to justify the 
use of separate latitude bands in the analyses. Whether the 
limits of each latitude band are optimum or not is undeter- 
mined, but it was felt that any further reduction in the 
width of each band would reduce the number of soundings used 
in each regression equation and thereby have a detrimental 
effect on a statistical analysis. The differences between 
the corresponding coefficients from the three- and four -day 
sample regressions were generally insignificant; however, 
this conclusion will be further tested by examination of the 
fractional explainted variance accounted for by the appro- 
priate regression equations when applied to the data- 
samples . 

In both Tables 2 and 3, the coefficients corresponding 
to channels 3, 4, and 5 appear to be weighted more heavily 
than any other triad of coefficients, which appears to be 
consistent with the response levels of Figure 1. 

A useful statistic for the fractional unexplained 
variance is 



where k is the number of predictors and n the sample 
size. 



2 



n - k - 1 
n - 1 




(9) 
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Thus the fractional explained variance is easily 
obtained from Equation (9) as v 

r 2 = 1 - S.-L... k . ~ 1 / std. e rr o r V (1 , 

n - 1 l std. dev. I 

y 

Values for R and of standard errors for each of the 
dependent tests corresponding to Tables 2 and 3 are listed 
in Table 4. 

TABLE 4. R 2 and Standard Errors (gpm) for NESC 
Dependent Regression Equations for 





500 mb and -300 
samples . 


mb 


Three- and Four 


-day 




500 


mb 






3 -Day Sample 


Band 1 




Band 2 


Band 3 


R 2 


.8050 




.8222 • 


.9050 


Std. error (m) 


61.72 




77.13 


69.68 


4-Day Sample 










R 2 


.8039 




.8154 


.9051 


Std. error (m) 


60.55 




76.78 „ 


69.79 




300 


mb 






3 -Day Sample 


Band 1 




Band 2 


Band 3 


R 2 


.87 90 




.. 8693 


.9351 


Std. error (m) 


77.51 




91.80, 


75.36 


4 -Day Sample 






‘V 

V • 


• 


R 2 


.87 64 




.8622 


.9349 


Std. error (m) 


76.22 




91.42 


75.50 
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As indicated in Table 4, Band 3 appears to produce higher 
values for R^, and therefore better results, than the other 
two bands in the NESC dependent regression analysis. The 
standard errors for the 300 mb tests were generally higher 
than the 500 mb tests, but one might expect this to happen 
since the standard deviations of Z N (3) are 9 reater than 
those of Z n ( 5) . In general, the 300 mb tests gave higher 
values of upon specification than the corresponding 

500 nib tests. 

B. DEVELOPMENT OF REGRESSION EQUATIONS USING 

FNWC- INTERPOLATED VALUES OF 500 MB AND 300 MB 
GEOPOTENTIAL HEIGHTS. 

Procedures similar to those used in Part A of this 
chapter were used on the FNWC 500 mb and 300 mb geopotential 
heights [Zp(5) and Zp(3) respectively] , corresponding 
to the identical SIRS scan-sounding positions and times. 

The same latitude band limits as specified in III. A. were 
used for both the three- and four -day regression analyses. 
The resulting coefficients for the FNWC regression equations 
are listed in Tables 5 and 6 for the 500 nb and 300 mb 
dependent tests, respectively. 
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TABLE 5. 500 mb FNWC Dependent Regression Coefficients 

for Three- and Four-day Data Base 







3 -DAY 


SAMPLE 


BASE 




Sample 


Size 


Band 1 
171 




Band 2 
190 


Band 3 
177 


C o 




5617.09 




5583.51 


4615.09 


C 1 




7.58 




-2.25 


-2.50 


c 2 




-19.82 




-2.58 


-9.43 


c 3 




23.50 




38.35 


57.19 


C 4 




25.55 




-6.14 


-16.06 


C 5 




-53.38 




-67.11 


-32. 17 


c 6 




-13.71 




5.85 


-9. 10 


c 7 




23. 14 




19.19 


12.22 


C 8 




-1.35 




-2.37 


0.80 






4 -DAY 


SAMPLE 


BASE 




Sample 


Size 


Band 1 
206 




Band 2 
227 


Band 3 
216 


Co 




5653.75 




5762.27 


4595.81 


Cl 




7.70 




0.36 


-2.24 


C 2 




-20.43 




-10.14 


-9.16 


C 3 




24.21 




46.38 


55.72 


C 4 




28.86 




-13.41 


-15.03 


c 5 




-59.44 




-64.52 


-30.50 


C 6 




-12.29 




4.74 


-9.58 


c 7 




27.24 




19.12 


8.71 


c 8 




-5.01 




-2. 15 


3.26 
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TABLE 6. 300 mb FNWC Dependent Regression Coefficients 

for Three- and Four-day Data Base 



3 -DAY SAMPLE BASE 



Sample Size 


Band 1 
171 


Band 2 
190 


Band 3 
177 


Co 


9156.82 


8451.37 


7337.29 


C 1 


14.03 


2.99 


* 


c 2 


-33.45 


-8.12 


-18.36 


C 3 


28.34 


34.19 


71.49 


C 4 


50.39 


29.29 


-3.58 


c 5 


-72.67' 


-86.77 


-35.32 


C 6 


-29.22 


5.55 


-9.36 


C 7 


21.93 


15.05 


9.18 


c 8 


8.82 


0.78 


* 




4 -DAY SAMPLE 


BASE 






Band 1 


Band 2 


Band 3 


Sample Size 


206 


227 


216 


c o 


9226.42 


8702.09 


7289.23 


C 1 


15.01 


7.09 


* 


C 2 


-36.75 


-20.98 


-16.37 


C 3 


31.79 


49.50 


66.94 


C 4 


52.29 


16.47 


* 


C 5 


-77.84 


-83.39 


-33.79 


C 6 


-27.24 


4.64 


-9.98 


C 7 


26.96 


15.72 


5.07 


C 8 


3.44 


0.41 


2.85 



*Indicates an insignificant contribution to the analysis of 
variance by this predictor. 
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The FNWC regression coefficients listed in Tables 5 and 
6 lend support to the observations from the NESC regression 
coefficients, in that the coefficients from the FNWC 
dependent test vary sufficiently between latitude bands to 
justify using separate bands. Furthermore, the differences 
between the three-day samples and four-day samples are again 
small; and the coefficients corresponding to Channels 3, 4, 
and 5 appear to be weighted more heavily than the others. 

Table 7 lists the values for r 2 and standard errors 
for each of the FNWC dependent tests listed in Tables 5 and 
6. As in the NESC tests, the values for R were higher in 
Band 3, standard errors were generally larger for the 300 mb 
tests than for the 500 mb tests, and the values for R^ 
were higher for the 300 mb tests than for the 500 mb tests. 

A comparison of the three- and four -day sample regressions 
for both the 300 mb and 500 mb cases revealed no significant 
differences in corresponding values of r 2, when subjected to 
to a t -test. 
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TABLE 7. 



R 2 and Standard Errors (gpm) for FNWC 
Dependent Regression Equations for 
500 mb and 300 mb Three- and Four -day 
Dependent Samples. 

500 mb 



3 -Day Sample 


Band 1 


Band 2 


Band 3 


R 2 


.7910 


.7954 


.9014 


Std. error (m) 


62.53 


84.55 


72.48 


4 -Day Sample 


R 2 


.7787 


.7899 


.8989 


Std. error (m) 


63.79 

300 mb 


83.69 


73.21 


3 -Day Sample 


Band 1 


Band 2 


Band 3 


R 2 


.8575 


.8361 


.9235 


Std. error (m) 


84.38 


105.99 


84.32 


4 -Day Sample 


R 2 


.8459 


.8328 


. 9222 


Std. error (m) 


85.77 


104.26 


85.07 
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C. COMPARISONS BETWEEN NESC AND FNWC DEPENDENT REGRESSION 
ANALYSES 

If one compares the results from Tables 2, 3, 5, and 6, 
channels 3, 4, and 5 appear to be weighted more heavily than 
the other channels for both the NESC and FNWC dependent 
tests. Table 8 combines the results of Tables 4 and 7 to 
facilitate comparison of NESC and FNWC dependent regression 
test results. Here, only comparative values of R^ by the 
various dependent stratifications are considered. If 
standard errors are required, reference should be made to 
Tables 4 and 7. 

Band 3 gave the highest values of R , averaging 8.7% 
higher over all dependent regression systems than Bands 1 
and 2, varying from 5.6% to 12.0% higher for individual 
comparisons. The 300 mb values of R^ averaged 4.65% 
higher overall than the 500 mb values of R^. with the NESC 
300 mb values of r 2 averaging 5.01% higher than the NESC 
500 mb values, and the FNWC 300 mb values of R^ averaging 
4.38% higher than the FNWC 500 mb values. 

Although quite close in most instances, in all cases the 
NESC values of R^ were slightly higher than the correspond- 
ing values of R^ from the FNWC tests, averaging 2.01% 
higher for the NESC tests. A closer look reveals an average 
of 1.69% smaller values, or "shrinkage," for the FNWC 500 mb 
tests as compared to the corresponding NESC tests. An 
average shrinkage of 2.32% occurred for the 300 mb FNWC 
'tests as compared to the NESC tests. The individual 
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shrinkages of for each band, averaged over both the three- 
and four-day data-samples is shown for the FNWC tests rela- 
tive to the NESC corresponding tests: 







Band 1 


Band 2 


Band 3 


500 


mb 


1.96% 


2.61% 


0.49% 


300 


mb 


2. 60% 


3.13% 


1.21% 


Bands 3 shows the 


smallest shrinkage 


of R 2 


values for FNWC as 



compared to the NESC dependent tests. That the NESC values 
of are higher than the FNWC values is not too surprising 
since the NESC data was originally the result of essentially 
an eight-predictor equation, the precise form of which was 
not known in this study. 

A comparison of three- and four -day overall regression 
systems showed that the three-day values of R were .46% 
higher than the four-day values. Since this is such a small 
difference, it is likely to be caused only by the additional 
data samples used in converting the particular three-day 
sample to a four -day sample, and the difference was not 
considered significant. 

Thus far, it is seen from the dependent-data tests on 
NESC and FNWC samples that the results are generally good, 
considering the r 2 values averaged .8579, with averages of 
.8679 for NESC and .8478 for FNWC. From the paragraphs 
preceding this one, it appears that the dependent tests 
produced better results in Band 3 than the other bands. 
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better results at 300 mb than at 500 mb and better results 
for NESC than for FNWC specification. 

It should be pointed out that the NESC and FNWC height 
data were highly correlated. The 500 mb values of R^ ranged 
from .9267 to .9699 and the standard errors varied from 35.26 
meters to 43.15 meters. The values of for the 300 mb case 
ranged from .9244 to .9707 and the standard errors varied 
from 51.46 meters to 59.64 meters. 
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IV. INDEPENDENT SAMPLE TESTS 



A. VERIFICATION PROCEDURE FOR THE INDEPENDENT DATA TESTS 

Each of the dependent regression equations, stratified by 
NESC and FNWC, by latitude bands, and for the different sizes 
of data-bases in Chapter III, was tested using independent 
data samples. The equations developed from the three-day 
samples were tested with independent data from a fourth day, 
for times centered at 12 hours and 24 hours later. The 
regression equations developed from the four -day samples were 
likewise tested for times centered 12 and 24 hours later 
using the equations developed from data collected over 24-27 
December, 1970. The corresponding four -day independent tests 
utilized data from the 28th (28/00 GMT and 28/12 GMT, i.e., 

12 and 24 hours later) . In addition, the 12- and 24-hour 
tests provided homogeneous statistics and were pooled to 
provide a larger number of independent data samples for each 
band and height. 

The results of the independent tests generate equations 
of the general form 





Z n (5,3) 


- D Q + D^Z n (5 , 3) 


(ID 


and 


Zf (5,3) 


= F 0 + F 1 Z p (5,3) 


(12) 



where Z N (5,3) and Zp(5,3) represent the height-estimators 
from the regression equation (Equation 8) , D Q and F Q are 

regression constants, and D-j_ and F-, the regression 
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coefficients. Zp(5,3) are the observed independent-test 
heights determined by interpolation from the data tapes of 
FNWC for the 500 mb and 300 mb levels. All independent 
Zn( 5,3) were the result of application of the NESC update 
regression to the independent set of scan-spot corrected 
radiances. Regression analyses for these tests were similar 
to the regression techniques applied in Chapter III, in that 
the newly computed height Z is now being compared by 
regression methods to the known height at each scan-sounding 
location. The independent" test results are presented in a 
slightly different manner in Tables 9, 10, 11, and 12. The 
mean and standard deviation for the NESC and FNWC verifica- 
tions, and the values for and standard error resulting 
from Zjj( 5) , Z N (3), Zp(5), and Zp(3) are listed in Tables 9, 
10, 11, and 12 for the three- and four -day independent tests. 
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TABLE 9 



. Verification Resulting From the Use of 

Equation (8) Upon 500 mb NESC Independent 
Data. 



3 -DAY TEST 









Sample 

Size 


Mean 


Std. 

Dev. 


R 2 


Std. 

Error 


Band 


1 


12 hr 


23 


5753.26 


99.04 


.7983 


45.52 






24 hr 


35 


5738. 60 


121.72 


.7 867 


57.07 






12+24 hr 


58 


5744.41 


' 112.59 


.7852 


52.64 


Band 


2 


12 hr 


31 


5464.45 


100.76 


.5428 


69.30 






24 hr 


37 


5477.38 


154.80 


.7652 


76.08 






12+24 hr 


68 


5471.48 


132.15 


.6934 


73.62 


Band 


3 


12 hr 


32 


5283.97 


215.61 


.8781 


76.54 






24 hr 


39 


5293.51 


227.12 


.9015 


72.24 






12+24 hr 


71 


5289.21 


220.49 


.8859 


75.01 










4 -DAY TEST 














Sample 

Size 


Mean 


Std. 

Dev. 


R 2 


Std. 

Error 


Band 


1 


12 hr 


28 


5760.27 


90.95 


.7537 


46.00 






24 hr 


57 


5739.58 


121.90 


.8308 


50.59 






12+24 hr 


85 


5746. 36 


112.52 


.8151 


48.67 


Band 


2 


12 hr 


33 


5513.99 


140. 14 


.7093 


76.77 






24 hr 


50 


5432. 87 


164.26 


.7499 


83.00 






12+24 hr 


83 


5465.11 


159.32 


.7509 


80.00 


Band 


3 


12 hr 


26 


5290.45 


202.12 


.7685 


99.26 






24 hr 


47 


5286. 16 


172.66 


.8281 


72.38 






12+24 hr 


73 


5287.69 


182. 30 


.7959 


82.94 
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TABLE 10 



Verification Resulting from the Use of 

Equation (8) Upon 300 mb NESC Independent 
Data. 



3 -DAY TEST 









Sample 

Size 


Mean 


Std. 

Dev. 


R 2 


Std. 

Error 


Band 


1 


12 hr 


23 


9427.13 


166.40 


.8980 


54.39 






24 hr 


35 


9398.46 


190.90 


.8526 


74.40 






12+24 hr 


58 


9409.82 


180.62 


. 8658 


66.76 


Band 


2 


12 hr 


31 * 


8941. 29 


157.58 


.7172 


85.23 






24 hr 


37 


8981.95 


202.08 


.8027 


91.04 






12+24 hr 


68 


8963.41 


182.96 


.7704 


88.34 


Band 


3 


12 hr 


32 


8682.44 


293.24 


.9084 


90.23 






24 hr 


39 


8695.46 


303.44 


.9334 


79.33 






12+24 hr 


71 


8689.59 


296.83 


.9207 


84.19 








4 -DAY TEST 














Sample 

Size 


Mean 


Std. 

Dev. 


R 2 


Std. 

Error 


Band 


1 


12 hr 


28 


9418.45 


160.38 


.8348 


66.42 






24 hr 


57 


9409.62 


203.85 


.8718 


7 3.64 






12+24 hr 


85 


9412.46 


189.70 


.8615 


71.03 


Band 


2 


12 hr 


33 


9012. 87 


206.88 


.7994 


94.15 






24 hr 


50 


8923.04 


234.25 


.8137 


102.15 






12+24 hr 


83 


8958. 68 


226.81 


.8148 


98.22 


Band 


3 


12 hr 


26 


8696.30 


266.23 


.8301 


111.99 






24 hr 


47 


8696.93 


242.90 


.8827 


84.13 






12+24 hr 


73 


8696.62 


249.60 


.8603 


93.94 
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TABLE 11 



. Verification Resulting From the Use of 

Equation (8) Upon 500 FNWC Independent 
Data. 



3 -DAY TEST 









Sample 

Size 


Mean 


Std . 
Dev. 


R 2 


Std. 

Error 


Band 


1 


12 hr 


23 


5736.63 


91.62 


.7572 


46.21 






24 hr 


35 


5724.33 


130.80 


.7103 


71.46 






12+24 hr 


58 


5729. 20 


116.11 


.7130 


62.75 


Band 


2 


12 hr 


31 


5468. 17 


123.09 


.6163 


77.56 






24 hr 


37 


5473.50 


159.20 


.7401 


82.31 






12+24 hr 


68 


5471.07 


142.86 


.6950 


79.49 


Band 


3 


12 hr 


32 


5275.86 


208.00 


.8440 


83.52 






24 hr 


39 


5291.72 


229.58 


.8863 


78.45 






12+24 hr 


71 


5284.57 


218.72 


.8627 • 


81.62 








4- 


-DAY TEST 














Sample 

Size 


Mean 


Std. 

Dev. 


R 2 


Std. 

Error 


Band 


1 


12 hr 


28 


57 50.87 


80.03 


.8329 


33.34 






24 hr 


57 


5728. 56 


120.04 


.8054 


53.43 






12+24 hr 


85 


5735.90 


108.52 


.8104 


47.53 


Band 


2 


12 hr 


33 


5522.73 


144.47 


.67 30 


83.94 






24 hr 


50 


5427.71 


171. 16 


.7143 


92.43 






12+24 hr 


83 


5465.49 


166.85 


.7170 


89.31 


Band 


3 


12 hr 


26 


5280.39 


198.52 


.6612 


117.93 






24 hr 


47 


5274.83 


180.50 


.7782 


85.95 






12+24 hr 


73 


5276.81 


185.76 


.7297 


97.26 
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TABLE 12 



. Verification Resulting From the Use of 

Equation (8) Upon 300 mb FNWC Independent 
Data. 



3 -DAY TEST 









Sample 

Size 


Mean 


Std. 

Dev. 


R 2 


Std. 

Error 


Band 


1 


12 hr 


23 


9396.04 


160.42 


.8465 


64.34 






24 hr 


35 


9378.17 


195.55 


.7 643 


96.36 






12+24 hr 


58 


9385.25 . 


181. 16 


.7883 


84.10 


Band 


2 


12 hr 


31 . 


8934.96 


187.09 


.7281 


99.23 






24 hr 


37 


8975.52 


217.55 


.7903 


101.03 






12+24 hr 


68 


8957.03 


203.76 


.7673 


99.04 


Band 


3 


12 hr 


32 


8665. 11 


282.86 


.8915 


94.70 






24 hr 


39 


8694.04 


308.89 


.9157 


90.88 






12+24 hr 


71 


8681.00 


295.70 


.9022 


93.13 








4 


-DAY TEST 














Sample 

Size 


Mean 


Std. 

Dev. 


R 2 


Std. 

Error 


Band 


1 


12 hr 


28 


9410.61 


144.12 


.7820 


68.57 






24 hr 


57 


9392.98 


200.90 


. 8310 


83.34 






12+24 hr 


85 


9398.79 


183.45 


.8198 


78.34 


Band 


2 


12 hr 


33 


9005.81 


217.27 


.7295 


114.81 






24 hr 


50 


8913.67 


244.64 


.7718 


118.07 






12+24 hr 


83 


8950.30 


237.16 


.7641 


115.91 


Band 


3 


12 hr 


26 


8678. 84 


267.39 


.7858 


126.31 






24 hr 


47 


8679.63 


257.69 


.8367 


104.31 






12+24 hr 


73 


8679.34 


259.33 


• .8188 


111. 18 
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B. SUMMARY OF INDEPENDENT TEST RESULTS 



1. Results Internal to the Independent Test 

Table 13 lists the combined 12 and 24 hour 
independent-test values of R^ abstracted from Tables 9, 

10, 11, and 12 to facilitate comparisons of test-methods 
contained in those tables. Again our emphasis is on R^ as 
a measure of the efficiency of the verification since this 
statistic has the effect of normalizing the variable standard 
deviations in Tables 9, 10, 11, and 12. 

Band 3, as in the dependent tests, gave the highest 

values of R^, averaging 6.99% higher than the other bands. 

The 300 mb values of R^ , again, as in the dependent cases, 

were higher than the 500 mb values, averaging 5.83% higher. 

The NESC values of R^ , as compared to the FNWC values of R^ , 

were higher, as before, by an average of 3.60% over all of 

the independent tests. However, the amount of shrinkage 

between the NESC and FNWC within the independent tests was a 

minimum in Band 2, as compared to the minimum values of 

shrinkage in Band 3 for dependent tests. This phenomenom is 

due to the fact that the independent test of R N ^ fell off 

considerably in Band 2, whereas the corresponding shrinkage 

was somewhat less with FNWC independent tests in Band 2. The 

2 

three-day estimator tests produced average values of R 
that were .76% higher than the corresponding four-day values 
for R . This, again, as in the dependent tests, is a small 
difference and is considered to be due to the small sample- 
size, or to the particular data-sample added to obtain the 
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TABLE 13. Summary of NESC and FNWC Independent Regression 
Test Results Using Equation (8) . 
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500+300 



four-day sample, and therefore of little statistical 
significance . 

2. Comparison of Independent- With Dependent-Data 
Specification. 

From the independent tests, the same basic 
conclusions as for the dependent tests are to be drawn since 
the results are quite similar. The independent-test values 
of R 2 were generally good with an average of .8004, an 
average of .8184 for R 2 from the independent NESC tests 
and an average of .7824 for the FNWC tests. The highest 
values of r 2 were in Band 3, highest at 300 mb, and 
highest for the NESC tests. 

In proceeding from the dependent tests to the 
independent tests, some loss of accuracy is to be expected. 
Values of R 2 again will be used for comparisons. Table 
14 is derived from Tables 8 and 13 to present an easier 
comparison of the corresponding values of R 2 from the 
dependent tests and the independent tests, with the value 
of shrinkage indicated also. In most cases the shrinkage 
of the R 2 values from the FNWC dependent tests to the 
independent tests was slightly larger than the shrinkage 
for the similar NESC tests. The FNWC shrinkage averaged 
1.59% more than the NESC shrinkage. 

The first group of numbers in Table 14 are averaged 
over all bands and over both three- and four-day coeffi- 
cient estimators. In . this group, the values of R 2 for 
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the NESC 300 mb tests show the lowest relative shrinkage of 
the four shown in the group. This might have been expected 
since the 300 mb NESC test results from both the dependent 
and independent tests showed better results than the others. 

The second group of numbers in Table 14 is averaged 
over the 500 mb and 300 mb level independent tests, and 
over three- and four-day coefficient estimators Z. This 
second group shows relatively small shrinkages in Band 1. 

This could possibly be due to the fact that in low latitudes, 
the radiance correction p'rocedure is somewhat more consis- 
tent than in mid-latitudes, and may be affecting the amount 
of shrinkage in a favorable manner. 

The third group of numbers in Table 14 is averaged 
over all three latitude bands, and shows the relative 
shrinkages by use of the Equation (8) based upon the three- 
and four -day dependent coefficient matrices. The amount of 
shrinkage of the three-day independent tests over the four- 
day tests, or vice-versa, is not very obvious, even though 
one might expect that the four -day tests should be more 
representative due to a larger input of data. From the 
findings of this study, it has consistently been determined 
that tests on independent data using the four-day data-base 
estimators are on the average slightly less effective in 
R.2 than when the corresponding three -day data base was 
used. 
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TABLE 14. 



Summary of Dependent to Independent 
Test Shrinkage. 



Test Being 
Compared 



Dependent 
Value of R 2 



Independent 
Value of R 2 



Average 

Shrinkage 



z n< 5 > 


.8428 


.7878 


5.50% 


Z F (5) 


.8259 


.7547 


7.12% 


Z N (3) 


.8929 


.8489 


4.40% 


Zp (3) 


.8697 


.8101 


5.96% 
avg . 5.7 5% 


Zjj Band 1 


.8411 * 


.8320 


.91% 


Zp Band 1 


.8183 


.7829 


3.54% 


Zjj Band 2 


.8423 


.7572 


8.51% 


Zp Band 2 


.8136 


.7359 


7.77% 


Z N Band 3 


.9201 


.8657 


5 . 44% 


Zp Band 3 


..9116 


.8283 


8.33% 
avg . 5.75% 


Zn( 5) 3-day 


.8441 


.7 882 


5.59% 


Z n ( 5) 4-day 


.8415 


.7873 


5.42% 


Z n ( 3) 3 -day 


.8945 


.8523 


4.22% 


Z n ( 3) 4-day 


.8912 


.8455 


4.57% 


Zp (5) 3-day 


.8293 


.7569 


7 . 24% 


Zp (5) 4-day 


.8225 


.7524 


7.01% 


Zp (3) 3-day 


.87 24 


.8193 


5.31% 


Zp(3) 4— day 


.8670 


.8009 


6.61% 
avg. 5.75% 
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V. CONCLUSIONS 



It has been shown in this study that regression equations 
developed from relatively small data bases of heights and 
corrected SIRS-B radiances can be used to compute heights on 
independent samples with relatively good results. It is felt 
that further testing should be conducted at other levels 
using comparisons of other samples to test the results 
obtained here. If the resources are available, larger data 
bases should be tested in a similar manner. Direct applica- 
tion of similar procedures could be of possible use in the 
update of FNWC analyses for standard levels in areas of 
sparse data, or even aboard ships with weather offices to 
provide real-time height calculations over an operational 
area. In the latter case the advantage of a small data base 
lies in the relatively small amount of core storage necessary 
for the type of computer calculations possible in shipboard 
operations. 

Although the NESC height data provided slightly better 
regression results than the FNWC height data, this was to be 
expected since the NESC heights were computed using an eight- 
predictor regression equation, and therefore were internally 
more consistent. The results obtained by using FNWC data 
for dependent and independent regression analyses, parallel- 
ing the procedures applied to the NESC data, lends support 
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to the overall concept of developing regression equations to 
compute geopotential heights using SIRS-B corrected radiances. 

Improvements of the results should be expected if 
longer-term data samples are used. Improvements should also 
be expected if the width of the latitude bands were smaller, 
and if the radiance data were subjected to gross error checks 
before entering into regression analyses. Also, a fact to 
be considered is the tendency of the SIRS-B channel response 
to degrade and require updated calibration after being in 
orbit for a period of time'. If degradation errors were 
linear in time, the regression-climatology should provide 
update corrections, but the normal case is that degradation 
occurs non-uniformly in the different channels. The best 
scan-soundings are normally obtained soon after the satellite 
is launched into its orbit about the earth, before appreci- 
able system-degradation can occur. 
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APPENDIX 



GRIDDING AND INTERPOLATION SCHEME 
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